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Abstract 

A distance measure is presented between two unitary propagators of quantum systems of differing 
dimensions along with a corresponding method of computation. A typical application is to compare 
the propagator of the actual (real) process with the propagator of the desired (ideal) process; the 
former being of a higher dimension then the latter. The proposed measure has the advantage of 
dealing with possibly correlated inputs, but at the expense of working on the whole space and not 
just the information bearing part as is usually the case, i.e., no partial trace operation is explicitly 
involved. It is also shown that the distance measure and an average measure of channel fidelity both 
depend on the size of the same matrix: as the matrix size increases, distance decreases and fidelity 
increases. 

1 Introduction 

Much of engineering design is predicated on having a model of the process which captures the essential 
features {e.g., the design variables and sources of uncertainty), a model of the desired process, and a 
means of comparing the two. In addition, it is also desirable that the same means of comparison can 
be determined or validated from experimental data. Distance and fidelity measures are proposed in 
[GLN05] for comparing what are referred to as real and ideal quantum processes represented as quantum 
operations [NCOO]. Several measures ai^e proposed in [GLN05] along with a set of criteria which aim at 
making the measure tractable both theoretically as well as experimentally. Here we focus more on the 
theoretical aspect and in particular using unitary propagators corresponding to the real and ideal systems. 
The proposed distance measure does, however, have an input output interpretation, and thus, in principal, 
it may be amenable to an experimental calculation although that is not pursued here. 



2 Problem formulation 

Consider a finite dimensional, closed, bipartite quantum system consisting of a part, S, referred to as the 
system, which bears the quantum information of interest, coupled to a part B, representing the bath or 
environment. The Hilbert space of the bipartite system is Hsb = 'Hs®'Hb where Tis and are finite 
dimensional with ns = dim'H^, ub = divnT-LB, and hence, nsns = dmiT-LsB- Let U G Q^snB^nsnB 
denote the unitary propagator of the (real) bipartite system, and let G G C"sx"s denote the unitary 
propagator of an ideal system acting only on Hs- The question addressed here is this: 

How close is U to G? 

Since it is implicitly assumed that experimental access is confined to the S'-system, it follows that per- 
formance is judged only by how well the reduced state output mimics the desired behavior in Hs as 
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represented by G. If, for example, U decomposes into a tensor product of unitaries, i.e., U = G ® Ub, 
then the S and B systems do not interact and the system is indistinguishable from G acting on Hs. This 
motivates expressing U as,^ 

U = G(^^ + R (1) 

where $ G Q'"-B>^nB |g unitary. Since any choice of $ determines R, it follows that it is always possible 
to express U in this form. Hence, consider the following distance measure between U G C"^'*^^"^'*^ 
andGG C"s^"s. 

d(C/,G) =imn| ||C/-G(g)$|| = Is } (2) 

where || || is any matrix norm on C"^"^^"'^"^. The term inside the norm is R as defined in (1). 
The unitary $ G C^sxriB jg yggd solely to find the closest propagator to U in the subspace 



G^^ 



$t$ = I where "closeness" is determined by the choice of norm. 
It is clear that using (2) as a distance measure does not consider the partial trace operation which 
is implicit in other distance (and fidelity) measures via the operator-sum-representation (OSR), e.g., 
[NCOO], [KW04], [GLN05]. As a result (2) is in general a conservative measure. More specifically, (2) 
is measuring distance over the whole space Hs x Ti-B rather than just over Hs- For example, suppose 
the (possibly correlated) pure state is acted upon by U and its output is compared to the ideal output 
of G (g) $ acting on \ip). The output state error is, 

e = (?7-G®#)|V') (3) 

It follows that 

max e'e = \\U - G 'S>MI (4) 

{i,\i,)=i 

where || -112 is the induced two-norm, i.e., the maximum singular value of the matrix argument. Suppose 
the input is random in the sense that |^) = with probabiUty pi for z = !,...,£. Let E denote the 
expected value operator taken with respect to the input distribution. Then, 

max E e'e = max Tr (U - G <^ ^)Uu - G p 

max Tr (U - G (gi ^)Uu - G (gi ^) p (5) 

p>0, Tr p=l 

= \\U - G <P\\l /nsnB 

1 /2 

where ||- Hp is the Frobenius matrix norm: ||X||p = {TrX^X)^^^ = (^^-j \Xij\^^ = (Ei 
with ai{X) the z-th singular value of X. Since the error is defined over states operating on the whole of 
HsxHb, (4)-(5) are upper bounds on the error when the actual input is uncorrelated, that is, when either 
IV') = 1^5) ® \iPb) or p = ps pb- Despite this, however, there is a benefit: using either the maximum 
singular value or the Frobenius norm in (2) gives an indication of both the effect of environmental input 
uncertainty (including unwanted correlations with the information input) and the "distance" of the actual 
system to an ideal. In the remainder of this note we will show how to compute the distance measure 
using both of the above matrix norms. We also explore the relation of this distance measure to a typical 
fidelity measure. 

'if the S and B channels are ordered reversely, then U = ^ <S) G + R. 
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3 Computing the distance measure with the Frobenius norm 



Consider the distance measure 

dF{U, G) = min 



\U-G®^\ 



(6) 



y/2nsnB 

where 1 / ^/2nsnB is a normahzation factor whick keeps dp (f/, G) in the range [0,1]. As shown below, 



dF{U,G) 



nsriB 



iTr 



1/2 



(7) 



G C 



where ||-||rp|. is the matrix trace-norm,^ the {C/(y)} are x matrix partitions of U, 



U = 



U, 



(11) 



(nsl) 



^(Ins) 



g QUsnB^nsUB 



(8) 



and {Gij} are the (scalar) elements of G, 

Gil 



G 



G 



Ins 



G. 



G 



nsns 



(9) 



Proof 

From the definition (6), 
1 



2nsnB 



U-G®^\ 



= 1 - 



nsriB 
1 

nsUB 
1 

nsns 



Re Tr (G^ ® $^)J7 



ij ^(ij) 



ns 

ReTr ^ G*,U, 
jj=i 



Re Tr T^^ 



(10) 



The first line above follows from the fact that [/, G, $ are all unitary. The last two lines 
follow from the definition of the tensor product together with (7)-(9). To compute dF{U, G) 
is equivalent to finding the maximum value of Re Tr over all unitary $ G C'*^^"^. A 
singular value decomposition of V gives V = WSV^ with unitary W and V and with the 
singular values in the diagonal matrix S = diag((Ji, . . . ,(t„^) G R"^, cti > (T2 > • • • > 
> 0. This gives, Re Tr r$t = Re Tr WSV^<^^ = Re Tr SV^^W = Re Tr 5.B 
with B = V^^^W. Clearly $ is unitary if and only if B is unitary. Hence, Re Tr SB = 
Y17=i ^ill achieve the maximum, ReTr5 = TrS' = ||r||rpj, when Re S = /, if 

and only if B = I, which is of course unitary. This completes the proof of (7). 



^For X e C" 



\x\\ 



= Tr VXtX = Yl 



min(n,m) 
i=l 



ai{X) where cri{X) is the i-th singular value of X. 
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Reversing the channel ordering 

By an analogous argument, if the S and B channels are ordered reversely from (1), then the distance 
measure becomes, 



dF{U,G) = min 



' in,.)" 



nsuB 

Tij = Tr G^'U^ij), i,j = l,...,nB 
where now { ?7(jj) } are x ns matrix partitions of U, 



U 



(11) 



u. 



(nsl) 



[riBriB) 



g QHsnsXrisnB 



Exact tensor product 

Suppose U is exactly a tensor product, i.e.. 



Then, 



U = Us®Ub, 

Unj) = (UshUB 

r = U,UGUUshUB={TrG^Us)UB 

|r||^ = TrJ\T^Gnj^fU^=\TrG^Us\TrlB = nB\TrG^Us\ 



The distance measure becomes, 

dF{Us0UB,G) 



1 

ns 



Tr G^Us 



1/2 



(11) 



(12) 



(13) 



(14) 



(15) 



If U and G have the same dimensions, that is, ne = 1, then the distance is zero if and only if U (= Us) 
and G differ by a scalar phase, i.e., Us = e^^G for some real phase angle cf). 



4 Computing the distance measure with the maximum singular value 

Consider the distance measure 



d2iU,G) = mm 



-L||c_o«*ll, 



(16) 



where l/\/2 is a normahzation factor which keeps d2{U, G) in the range [0, 1]. In this case we have not 
been able to obtain an exact solution, but we can estabhsh the following bounds: 



^ "U -G(^^L<d2{U,G) < ^ 



V2 



V2 



U-G0^ 



(17) 
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where $ is the solution to the convex optimization problem: 



minimize ||J7 — G (8) <&||2 
subject to < Ib 



(18) 



and where $ is obtained from $ via the singular value decomposition of $ as follows: 



(19) 



The upper and lower bounds in (17) will be close to each other if the singular values of $ are close to 
unity, that is, if $ is close to a unitary. 



Problem (18) is a convex optimization because any norm is a convex function, the argument 
in the norm is affine in the optimization variable and the constraint set < Ib is 
convex [BV04]. However, there is no guaranty that the resulting optimizer $ is unitary, i.e., 
on the boundary of the constraint. Since $ in (18) is less constrained then in (16), it follows 
that the lower bound in (17) applies. In consequence, consider $ from (19) as a unitary 
approximation to and since it is not necessarily the optimal solution to (16), the upper 
bound in (17) follows which completes the proof. 

5 Application to optimal control or system design 

In the ideal case where U and G are of the same dimension, the term |Tr G^Us\ in (15) has often been 

proposed as an objective function to be maximized for optimal control design (e.g., [PK02]) and also 
for determining properties of the control landscape [RHR04]. This is equivalent to posing the distance 
measure (15) as an objective function to be minimized. Where U and G are not of the same dimension, 
the distance measure (7) has been reported in [GBR+] for optimal control design. 

In general, for either control design or system design, U will depend on some parameters. In many 
cases these parameters are constants, e.g., coefficients of specified time functions which make up the 
control field, settings of wave-plate angles in a photonic device, geometry variables in a circuit layout, 
and so on. In these cases, U = U{9) where is a (constant) vector to be selected out of a set 9 to 
minimize the distance measure (2). Equivalently, consider the following optimization problem: 



where || || is any matrix norm and the optimization variables are $ and 9. Although typically 9 is a 
convex set, e.g., & = { ^ | ||^ — ^o|| < <^ }» in general (20) is not a convex optimization. In the first 
place, it is almost never the case that U{9) is a convex function of 9. Suppose for example, that the 
Hamiltonian which generates U{9) depends linearly on 9, i.e., H{9,t) = 9k Hk{t). If U{9) is the 
associated propagator which makes its appearance at some time r, then it is unUkely that U{9) is convex 
over 9. Secondly, the constraint = 7b is not a convex set. This constraint, however, can be relaxed 
to the convex set < Ib resulting in. 



Since \\U{9) — G(^^\\ for any norm is a convex function of $ and U{9), and < 7b is a convex 
set in it follows that (21) is a convex optimization over U{9) and $ (see, e.g., [BV04, §3.25]). Again, 



Proof 



minimize \\U{9)-G'^^\\ 
subject to = 7b, 6' G 9 



(20) 



minimize ||;7(6') - G O <^|| 
subject to < 7b, 9 e@ 



(21) 
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this is not a convex optimization over because U{0) is not a convex function of 9. Nevertheless, the 
following iterative scheme will always find a local solution by reducing the distance measure in every 
step. 



Initialize 
Repeat 

Until 



= 60 

1 . $ = arg min < 

2. 9 = arg min 



U(9) - G (g) $ 



\U{9)-G(g)^ 
U{9) - G (8) $ 
stops decreasing, 



In step 1 if the Frobenius norm is used then $ can be calculated exactly (7). Under the two-norm we 
would use the approximation from (19). Step 2 requires using a local solver. 



6 Distance and fidelity 

We will show that the distance measure defined here and a typical measure of fidehty both depend on 
the size of the matrix V as given by either (7) or (11). There are many ways that fidelity has been 
defined to compare the desired unitary G with a quantum channel, e.g., [NCOO], [KW04], [GLN05]. 
Specifically, let S denote a trace-preserving quantum channel mapping states ps € Hs to p5 G Hs with 
the operator-sum-representation (OSR), 

ps = S{ps) = Y,SkPsSl Y.SlSk = Is (22) 

k k 

with (matrix) operation elements € C"^^"^. Consider, for example, the worst-case pure state fidelity, 
/(<S,G) = min(G|V'5))^fe(G'|V5)) = minJ]|(V5|GtS,|V'5)|' (23) 

where = X^fe 5^ |V's')(V's I 'S'| is the reduced output state of 5 with l^/;^) the pure input state. In general 
calculating f{S, G) is not easy as it is not a convex optimization; it is, however, bounded as follows: 

/(<S,G)</(5,G)<7(5,G) (24) 

where 

k 

(25) 

f{S,G) = minV|TrGt5feP5|' 
— PS ~^ 

k 

All these fidelities are in [0,1] and equal to one if and only if S{ps) = GpsG^ for all densities ps- Both 
bounds are easy to obtain: f{S, G) is direct and f{S, G) is a convex optimization over all densities {i.e., 
over all ps G C"'^^"^, > 0, Tr p^ = 1), and hence, can be numerically obtained. If the worst-case 
density associated with f(S, G) is nearly rank one, then f{S, G) ^ f{S, G). 

In order to relate U to S assume that the input state to U is the pure uncorrelated state^ \tpB) <X) \i>s) 
where \iPb) has elements ipBi, ^ = 1, • • • , with IV'BiP = 1- It then follows that 

riB 

Sk = ^i^BiU(ki), k = l,...,nB (26) 

1=1 

'For convenience we use what we referred to previously as the reversed channel ordering. 



6 



where the i7(jfcj) are the ng x ng matrix partitions of U as given by (12). The upper bound in (24) is then, 

1 '^'^ f \ 1 

with Tki = Tr G^Ufj^i^ from (11). This upper bound on pure state fidelity thus depends on V weighted 
by the state of the environment |V'b)- Suppose the environment state is completely random, that is, 
E{|?/'b)(V'b|} = Is/nB- Then the average value of the pure state fidelity upper bound becomes. 



E{/([/,G)} = ^— ir/cil' = ^— lirill (28) 

ns^UB ns^UB 

The same result holds with F defined by (7). Thus the distance measure presented here and a typical 
channel fidelity both depend on the matrix F defined by either (7) or (11), depending on how channel 
ordering is ascribed. Increasing the size of F increases fideUty and decreases distance. 

We also mention that the inequaUty relation of distance and fideUty described in [GLN05] also holds 
here as well, namely, 

(1 - di,iU, Gjf < E {m G)} < 1 - dF{U, Gf (29) 

Proof 

We first show that the bounds follow from ||F||2 = maxj cTj < ns where crj is the i-th 
singular value of F. Then we will show that ||F||2 < ns holds. 

The upper bound holds if and only if ||F||p < ns ||F||rpj., or equivalently, af < ns Oi. 
Assuming that crj < ng, it immediately follows that of < ns J2i '^i which estabUshes 
the upper bound. 

The lower bound is equivalent to ||F||rpj. + (l/ng) ||F||p < 2^^7X6 ||F|| p. Assuming again 

1 1 1 1 2 

that IIFII2 < ns, and hence, ||F||p < ns ||F||rpj,, the lower bound will hold if ||F||rpj. < 
^/nB llrllp- This is a known inequality for any matrix and thus establishes the lower bound. 

To show that cjj < ns, observe that since the elements of F (11) are obtained via the trace 
operation, then by definition they are equivalently obtained from the sum of the diagonal 
elements of G^U^ijy Let U' denote the matrix whose ns x ns matrix partitions, U'^^-y are 
diagonal, with diagonal elements the diagonals of G^U(^ijy Thus, F^ = Tr G^^Uf^ij^ = 
Tr U'^-jy Now replace the trace operation by the vector products, that is, Fjj = wfU'wj 
where Wi is an nsnB x 1 vector with ns ones in elements 1 + (i — 1)715 to ins and zeros 
elsewhere. Hence we can write F = W'^U'W where W = [wi ■ ■ ■ Wng] is nsnB x n^. 
Observe that \\W\\2 = s/ns and since U' is constructed from elements of unitary matrices 
Gand[/,then||;7'||2 < 1. It therefore follows that ||F|| 2 = ||W^?7'VF||2 < ||t^||2 WW^ < 
ns- This completes the proof. 
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